function h=ellipse(ra,rb,az,x0,y0,C)


t = 0:0.1:2*pi+0.1;

x = x0 + ra*cos(t)*cos(az)-rb*sin(t)*sin(az);
y = y0 + ra*cos(t)*sin(az)+rb*sin(t)*cos(az);

plot(x,y,C, 'Linewidth', 2);
